Baryonic resonances close to the KN threshold: the case of S(1385) + in pp collisions 
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We present results of an exclusive measurement of the first excited state of the E hyperon, E(1385) + , pro- 
duced in p + p — y E + + K + + n at 3.5 GeV beam energy. The extracted data allow to study in detail the 
invariant mass distribution of the E(1385) + . The mass distribution is well described by a relativistic Breit- 
Wigner function with a maximum at m = 1383.2 ± 0.9 MeV/c 2 and a width of 40.2 ±2.1 MeV/c 2 . The 
exclusive production cross-section comes out to be 22.27 ± 0.89 ± l.5&^2°vo A 1 ^- Angular distributions of the 
E(1385) + in different refer ence frames are found to be compatible with the hypothesis that 33 % of E(1385) 
result from the decay of an intermediate A ++ resonance. 

PACS numbers: 25.75.Dw,25.75.-q 
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The excitation spectrum of hadrons reflects proper- 
ties of quantum chromo dynamics (QCD) in the non- 
perturbative sector. Many hadrons, mesons as well as 
baryons, can be identified as members of multiplets 



2 



within quark models. More advanced techniques, such 
as various effective theories which are directly related 
to QCD, however, generate dynamically certain res- 
onances as emerging from interactions among other 
(more fundamental) hadrons |Q]]. This dual picture mo- 
tivates the quest for understanding the very nature of 
resonances or hadrons in general. 

Considering the S = — 1 strange baryon sector, the 
£(1385) is the first excited state of the £ hyperon and 
has a spin of 3/2 h and isospin 1, in analogy to the 
A(1232) as first excited state of the nucleon. This res- 
onance is characterized by a short life time that trans- 
lates into a natural width of 35.8 ± 0.8MeV/c 2 Jl. 
The £(1385) itself is considered as a standard quark 
triplet but its vicinity to the A(1405) in the mass spec- 
trum correlates the study and understanding of the two 
resonances. Indeed, as pointed out in 0], within a chi- 
ral SU(3) lagrangian approach the anti-kaon spectral 
distribution is intimately related to the A, A(1405), 
A(1520), £ and £(1385). In particular, the A(1405) 
arises naturally from hadronic effective models as a 
mixture of a K~ N and 7r£ bound systems JH-Q] but 
the contribution of the two 'molecular' states and the 
dependency of the A(1405) formation upon the reac- 
tion is still far from being understood. Any analysis of 
the lowest lying hyperon resonances suffers from the 
overlapping mass distributions, from possible interfer- 
ence effects through their K — N coupling, and from 
the common decay into £ — it. In a first step of such 
an analysis we will determine the £(1385) + produc- 
tion characteristics in the A — it decay channel (this 
paper) and address A(1405) production in a forthcom- 
ing publication. 

The largest part of the present knowledge on the 
£(1385) hyperon has been gained by employing 
photo-production on a proton target M\ (for the cor- 
responding theoretical analysis cf. J9]), K~p col- 
lisions iflOl [Till , and pp collisions Q2|. The latter 
measurement was accomplished at a beam momen- 
tum of 6GeV/c corresponding to an excess energy, 
defined as the energy above the production threshold, 
of e = 830 MeV. The HADES collaboration has re- 
cently measured pp collisions at 3.5 GeV beam en- 
ergy. The £(1385) + is identified in the K + npn~ n + 
final state via the Att + decay of the £(1385)+ (BR= 
88 %), where the A is identified by the pir~ decay 
(BR= 63.9 %). The present paper reports on this mea- 
surement. 

Thinking of the exclusive £(1385) + -fT + n produc- 
tion in pp collisions the one-boson exchange model 
can be considered. In this framework several diagrams 
are conceivable: (i) tt + , p + exchange and £(1385) + 
production in a ir + (p + )p — £(1385) + AT + vertex with 
possible excitation of an intermediate A ++ decaying 
into £(1385) + A'+; (ii) K°W exchange and produc- 
tion in the pK 0( -*~> - £(1385)+ vertex. Another pos- 
sibility is related to internal meson conversion in a 
tt + K + K°* vertex and £(1385) + production in the 



K°*p — £(1385) + vertex. Pion exchange diagrams 
have been considered in the analysis lfl2tl . while in lfl3ll 
kaon exchange has been included as well. Particularly 
interesting is the role of the intermediate A ++ exci- 
tation, as discussed in (Till . The angular distribution 
of the £(1385) + extracted in this work provides in- 
formation which can serve to disentangle the different 
production mechanisms. 

The results presented here have also to been seen 
in the context of a forthcoming analysis of the p(3.5 
GeV)+Nb reaction, investigated with the same appa- 
ratus, which allows for studies of the medium modi- 
fications of the £(1385). In the nuclear medium, the 
£(1385) is predicted to suffer an attractive interaction 
encoded in the real part of the self-energy resulting in 
a 40 MeV "mass shift" and a broadening of the width 
by a factor of two at normal nuclear matter density 
1151 . For further discussions of in-medium modifica- 
tions of £(1385) see EMI- A first identification of 
the £(1385) in sub-threshold heavy ion collisions has 
been reported in lfl9ll . in which the statistical error of 
the extracted signal did not allow to draw conclusions 
about a possible broadening of the resonance s pect ral 
function. The measured yield of the £(1385) IT19I1 is 
well reproduced by a statistical model. 

Baryon resonances, such as the A(1520) and the 
£(1385), have also triggered the interest in the con- 
text of heavy ion reactions at higher energies I20I, l2lll . 
There, the measurement of the resonance abundance is 
thought to reflect the dynamics of their production and 
can be connected with a stage of the reaction prior to 
hadronisation l22tl . 

Following this line of reasoning, we present here an 
analysis of the spectral shape of the £(1385) + recon- 
structed in the exclusive reaction p(3.5 GeV) + p — >• 
£(1385) + +A' + +n. Our work is organized as follows. 
Sections II and III present the experimental set-up to- 
gether with the event and particle selection. Section IV 
deals with techniques developed to identify the differ- 
ent background sources, and in section V the extracted 
£(1385) + signal is discussed. Sections VI and VII 
deal with the differential cross section and the mod- 
eling of the acceptance corrections, respectively. We 
close with a summary and conclusions in section VIII. 



II. THE EXPERIMENT 

The experiment was performed with the High 
Acceptance Di-Electron Spectrometer (HADES) at the 
heavy-ion synchrotron SIS 18 at GSI Helmholtzzen- 
trum fur Schwerionenforschung in Darmstadt, Ger- 
many. HADES is a charged-particle detector con- 
sisting of a 6-coil toroidal magnet centered on the 
beam axis and six identical detection sections (with 
nearly complete azimuthal coverage) located between 
the coils and covering polar angles between 18° and 
85°. Each sector is equipped with a Ring-Imaging 
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Cherenkov (RICH) detector followed by Multi-wire 
Drift Chambers (MDCs), two in front of and two be- 
hind the magnetic field, as well as two scintillator ho- 
doscopes (TOF and TOFino). Lepton identification 
is provided mostly by the RICH and supplemented at 
low polar angles with Pre-Shower detectors, mounted 
at the back of the apparatus. Hadron identification is 
based on the time-of-flight and on energy-loss infor- 
mation from the scintillators and the MDC tracking de- 
tectors. In the following the TOF-TOFino-PreShower 
system is referred to as Multiplicity And Time-of- 
flight Array (META). A detailed description of the 
spectrometer can be found in 12311 . 

A proton beam of ~ 10 7 particles/s with 3.5 GcV 
kinetic energy was incident on a liquid hydrogen tar- 
get of 50 mm thickness corresponding to 0.7% inter- 
action length. The data readout was started by a first- 
level trigger (LVL1) requiring a charged-particle mul- 
tiplicity, MUL > 3, in the META system. A total of 
1 . 14 x 10 9 events was recorded under these experimen- 
tal conditions. 

A common start time for all particles measured in 
one event was reconstructed as described in detail in 
rf25ll . With this method the time-of-flight information 
could be determined for almost all of the measured 
events. 



in. PARTICLE IDENTIFICATION AND EVENT 
SELECTION 

The aim of this analysis is to reconstruct the 
S(1385) + signal in the following reaction and decay 
chain: 

p(3.5GeV) + v S(1385)+ + K + + n (1) 

A + 7T+ 
L^. p + 7T~ . 

The value of BR= w 56 % accounts for the total 
branching ratio via this decay chain. 

The first analysis step consists of selecting events 
containing four charged particles (p, 7r~, tt + , K + ). 
Particle identification is performed employing the en- 
ergy loss (dE/dx) of protons and pions in the MDCs 
and adding the dE/dx information from the TOFino 
system to further select kaons. The kaon selection is 
performed by TOFino only, because only a negligi- 
ble amount of kaons are observed in the TOF region. 
Fig. Q] shows the dE/dx distribution as a function 
of the particle momentum extracted from the MDCs 
(panel (a)) and TOFino (panel (b)). The full black 
curves refer to adapted Bethe-Bloch formula 12411 and 
the dotted and dashed lines show the cuts applied 
to select pion and proton candidates in the MDC: 
(panel (a)) and kaons in the TOFino (panel (b)). The 



masses of these kaon candidates are calculated us- 
ing momentum and velocity. The purity of the kaon 
signal is then enhanced by selecting masses between 
280 - 780MeV/c 2 . The reconstruction of the A hy- 
peron is achieved exploiting the decay A — > p + tt~ 
(BR= 63.9 %). Since the A hyperon has a mean decay 
length of cr « 7.89 cm, selections on the decay vertex 
can be applied to reduce the background. 

The primary reaction vertex has been estimated by 
computing the point of closest approach between the 
7r + and the K + track pairs. The following topological 
cuts are employed to select A candidates and to en- 
hance the signal-to-background ratio in the p-ir invari- 
ant mass distribution: (1) distance between the pro- 
ton and pion tracks (dp-^- < 20 mm), (2) distance 
of closest approach to the primary vertex for the pro- 
ton (DCAp > mm) and for the pion (DCA W - > 
DC A p ), (3) distance between the primary reaction ver- 
tex and the secondary decay vertex (d(A — V) > 
15 mm). The resulting invariant-mass distribution of 




Momentum [MeV/c] 

FIG. 1: Color online. dE/dx versus momentum for all de- 
tected particles in the MDC (panel (a)) and TOFino (panel 
(b)) detectors. The solid black curves indicate adapted 
Bethe-Bloch formulas and the dashed ones show the inter- 
val used to select protons and pions in the MDC (panel (a)) 
and kaons in TOFino (panel (b)). 

the selected p-7r~ pairs is shown in Fig. [2 in which 
the contribution of the A hyperon is clearly visible. 
The background was fitted using a Landau function 
plus a polynomial of third order, while a Gaussian fit 
was applied to the signal. The obtained mean value 
and variance for the reconstructed mass are Mj\ = 
1115.2 McV/c 2 and a A = 2.5McV/c 2 , respectively. 
In order to further improve the data sample, only those 
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2: Color online. Invariant mass distribution of selected 
p-ir~ pairs after the K + selection (see text for details). The 
two vertical lines define the mass interval select events with 
A candidates. 



events are selected for which the p-n~ invariant mass 
is within 1110 - 1120 MeV/c 2 (corresponding to 2.0 
a), as indicated in Fig.[2]by the vertical lines. A further 
step consists in the selection of the missing neutron 
appearing in reaction ([T), which cannot be detected in 
HADES. The missing mass to the four charged parti- 
cles (p, ir~, 7r + , K + ) is calculated after the A and K + 
selection and the result is shown in Fig. [3] A clear 
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FIG. 3: Color online. Missing mass to the four charged par- 
ticles p, 7T~, 7r + and K + , after the selection of a A candidate 
in the same event. The gray full curve shows the sum of all 
identified contributions. 



signal close to the nominal neutron mass is visible 
(M„ = 939.5 MeV/c 2 , <r n = 12.6 Me V/c 2 ) stick- 
ing out of a broad background distribution. Simulation 
studies have shown that the background to the left of 
the neutron peak is entirely caused by events with pro- 
tons and pions misidentified as kaons. The red and the 



blue histograms in Fig. [3] show the background due to 
pion and proton misidentification obtained by a side- 
band analysis of the experimental data, as discussed 
in section [IV] Additionally, Fig. [3] illustrates the con- 
tribution to the neutron signal from the non-resonant 
channel p + p — >• A + K + + n + ir + (black histogram) 
and p + p — >• S° + K + + n + tt + (green histogram). 
While the shape of these two contributions is obtained 
by full-scale simulations, the absolute yield of these 
channels can be determined by fitting the total experi- 
mental spectrum. 

Fig. [4] shows the K + mass distribution after 
the A selection and a cut around the neutron 
mass 877MeV/c 2 < MM(p, K + , ir+, tt") < 
999 MeV/c 2 . This cut is shown by the vertical dashed 
lines in Fig.[3]and corresponds to ± 5 a. 
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FIG. 4: Mass distribution of the selected K + after the cuts on 
the A and missing neutron candidates (see text for details). 
The vertical lines represent the nominal masses of ir + , K + 
and proton respectively. The gray-shaded area shows the ap- 
plied selection on the K + mass. 



IV. BACKGROUND IDENTIFICATION 

Due to the incomplete purity of the kaon selection, 
even after all the cuts applied so far, a large fraction 
of the kaon candidates are still pions or protons. In 
particular, the contamination due to misidentified pi- 
ons is rather high as visible in Fig. [4] In order to 
estimate quantitatively this background contribution a 
dedicated side-band analysis based on the mass distri- 
bution of the selected K + candidates has been devel- 
oped. The goal of this analysis is to produce an event 
sample, which does not contain any kaon but either a 
pion or a proton misidentified as a kaon. To emulate 
the background, events have been selected containing 
four charged particles: one proton, one 7r~, one ir + 
and a fourth positively charged particle for which no 
identification is required. The mass distribution of this 
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fourth particle is shown in Fig. [5] where one can see 
that most of the yield is located around the nominal 
pion and proton masses and that no kaon peak is visi- 
ble. These events represent a sample with a dominant 
contribution of misidentification background. A data 
sample with misidentified pions is explicitly extracted 
from the distribution in Fig. [5] by selecting the mass 
range of —0.25 — (GeV/c 2 ) 2 , while for protons the 
interval from 0.615 to 4 (GcV/c 2 ) 2 is chosen. These 
background samples are analyzed in the same way as 
described in section|IIl]with an exception made for the 
secondary vertex cuts of the A decay. These cuts ((1- 
3) in sectionHIIli have been omitted in order to enhance 
the statistics of the background sample. 

However, although the selection of the misidenti- 
fied pions and protons is done in a range rather distant 
from the nominal kaon mass, the two misidentification 
samples still contain a certain amount of real kaons. 
These real kaons are produced together with a A, due 




x10 b 



Mass 



1 2 

[(MeV/c 2 ) 2 ] 



FIG. 5: Count rate of all K candidates as a function of 
the squared mass. The two gray-shaded areas show the se- 
lected 7r + and the protons on the left and right side, respec- 
tively. The dashed verticals line indicate the nominal squared 
masses of tv + , K + and protons. 



to strangeness conservation. Since the selection on the 
7r~-proton pairs invariant mass around the A nominal 
mass is also applied to the background samples, the 
contribution by real kaons is enhanced. 

Our strategy to reduce the contribution of the real 
kaon signal to this background samples is to smear si- 
multaneously the momenta of the identified tt + and ir~ 
such that the A signal in the n~ -proton invariant mass 
distribution disappears. This way, theA selection in 
the p-TT invariant mass distribution does not enhance 
the fraction of events with a K + . 

A further issue to be addressed deals with the spe- 
cific kinematics of the particles selected via the cut 
around the nominal kaon mass (from here on consid- 
ered as fake kaons) and of the explicitely misidenti- 
fied pions and protons, employing the side band cut. 



The fake kaons correspond to particles with the 'cor- 
rect' kaon mass but which are not real kaons. They 
are contributing to the missing mass distribution shown 
in Fig. [3] mainly in the range away from the nominal 
neutron mass and are hence a source of background. 
Figure [6] compares the momentum distribution of the 
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FIG. 6: Color online. Momentum distribution of the fake 
kaons (bullets) together with the protons (triangles) and pi- 
ons (squares) explicitly misidentified as kaons. 



fake kaons (bullets) together with the momentum dis- 
tributions of the misidentified protons (triangles) and 
pions (squares). One can see that the shapes of these 
distributions are rather different. This difference can- 
not be neglected during the evaluation of background 
due to fake kaons, since we are performing an exclu- 
sive analysis with the total four-momentum being con- 
served. 

In order to reproduce the momentum distribution of the 
fake kaons the following function must be determined: 



7(b1) 



f K+ (\p\) 



an+f*+(\p\) + fp(\p[)' 



(2) 



where fx+(|p|), f^+dpl) and f p (|p|) are the momen- 
tum distributions of the fake kaons, misidentified pi- 
ons and protons, respectively. The factor a„.+ accounts 
for the unknown relative contribution of the misidenti- 
fied pions and protons to the total background. Given a 
value of the parameter a w +, each event corresponding 
to the pion 'side' of the mass distribution in Fig. [4] is 
weighted by the factor a 7l + and multiplied by 7(|p|), 
while each event from the proton 'side' is weighted 
with 7(|p|) only. The total missing mass is calculated 
for each weighted event and the obtained distributions 
are compared to the background underlying the neu- 
tron signal. Figure [3] shows the resulting missing mass 
generated by the two background samples (blue his- 
togram for protons and red histogram for pions). The 
parameter a n + is varied systematically until the sum 
of the two background distributions gives the best fit 
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to the region of the spectrum on the left side of the 
neutron peak, i. e. for < MAl(p,K + ,ir + ,Tr~) < 
870McV/c 2 . Together with the such derived back- 
ground the contribution of the signal channels p + p — > 
A + K+ + n + 7T+ and p + p -> E° + K+ + n + n+ 
are displayed in Fig. [3] (black and green histograms re- 
spectively). The gray curve in Fig. [3] shows the sum 
of all contributions, which reproduces the experimen- 
tal data very well. Our side-band method allows in- 
deed to reproduce the total background. Furthermore, 
the contamination by misidentified pions and protons 
to the signal can be evaluated quantitatively. After 
the cut around the nominal neutron mass in the to- 
tal missing mass distribution one obtains precisely the 
signal-to-background ratio and hence can estimate the 
contribution of the misidentification background to the 
S(1385)+ spectrum. 



able, Tq the resonance width, T the mass-dependent 
width, i*i(q) the Blatt-Weisskopf parameter and R = 
1/197.327 Me V -1 the centrifugal barrier parameter. 
Accordingly, q, qo, m, mo, T and To are in MeV units 
here. The Blatt-Weisskopf centrifugal-barrier parame- 
ters l27ll absorb possible divergences of the energy de- 
pendent width. In order to quantify the different contri- 
butions to the spectrum in Fig. [7] a fit is applied to the 
data in the mass range 1250 McV/c 2 < M(A— 7T + ) < 
1530 MeV/c 2 (vertical dashed lines in Fig . [7J . The to- 
tal fit function is composed of: 

• the fixed distribution of the misidentification 
background and of channel (i), 

• the phase-space distribution of reaction (ii) fil- 
tered through the geometrical acceptance and 
the detector efficiency scaled by one fit parame- 
ter and 



V. THE S(1385)+ SIGNAL 

In the final step events conforming to the A, ir + , K + 
and neutron hypothesis are used to calculate the A-7r + 
invariant mass. The resulting distribution is shown in 
Fig. [7] A clear peak at the E(1385) + pole mass is vis- 
ible on the top of some background. The contribution 
of the misidentification background (blue histogram) 
is shown together with the contributions of the direct 
reactions (i) p + p -+ E° + n + + K + + n (green his- 
togram) and (ii) p+p -+ A + 7T+ + K + + n (black his- 
togram). After the determination of the scaling factors 
a 7T + and 7(|p|) the form and yield of the misidentifica- 
tion background is fixed. The reactions (i) and (ii) rep- 
resent a source of background with the same hadrons 
in the final state as the reaction (1). The absolute yield 
of the channel (i) is determined by fitting the signal in 
the neutron mass spectrum shown in Fig. [31 as already 
mentioned, while the yield of channel (ii) can only be 
determined by fitting the total distribution shown in 
Fig. |7] We chose to fit the experimental S(1385) + sig- 
nal with a Breit-Wigner function folded with the effi- 
ciency and acceptance corrections thus fitting the data 
instead of correcting the extracted experimental signal. 
This method allows us also to evaluate the contribution 
of the background source (ii) more precisely. The used 
p-wave relativistic Breit-Wigner function 112611 reads as 
follows: 



Breit-Wigner oc 



mlTl 



ql (m 2 - m 2 ) 2 + m 2 r 2 ' 
,3 



(3) 



r = r 



mpq" 



Fi(q) = 



1 + (qoR)' 
1 + (qR) 2 



where q is the momentum of the decay products in 
the S(1385) + rest frame, qo the momentum that cor- 
responds to the pole mass mo, m the mass vari- 



• a corrected relativistic Breit-Wigner function 
with fit parameters for the mass, width and 
height. 
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FIG. 7: Color online. Invariant mass distribution of the se- 
lected A-7r + pairs. The circles show the experimental data, 
the black and green curves depict the contributions coming 
from the non-resonant production, the blue curve shows the 
misidentification background and the gray dashed curve rep- 
resents a Breit-Wigner distribution fitted to the data. The 
solid gray histogram represents the result of the final fit (see 
text for details). Dashed vertical lines depict the interval cho- 
sen for the fit. 



The correction of the Breit-Wigner function is nec- 
essary to account for the finite available phase space 
and for the effects of the geometrical acceptance of 
the spectrometer and the efficiency of the analysis cuts 
that modify the measured spectral distribution in com- 
parison to a pure Breit-Wigner function. The geo- 
metrical acceptance and reconstruction efficiencies for 
the S(1385) + have been estimated employing a Breit- 
Wigner distribution of the resonance as the input of a 
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full-scale simulation. The procedure applied to extract 
the correct acceptance correction for the E(1385) + 
signal is explained in detail in sections IVIBI These 
corrections also account for the finite resolution of the 
apparatus which cause a broadening of the £(1385) + 
resonance of about 3MeV. Since the fitted Breit- 
Wigner distribution is filtered through this correction, 
the extracted width and pole mass are already cor- 
rected also for the resolutions effects. In addition, the 
phase-space correction factor as a function of the in- 
variant mass must be taken into account. This fac- 
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FIG. 8: Color online. Invariant mass distribution of the se- 
lected A — 7r + pairs after the background subtraction. The 
black circles show the experimental data together with the 
statistical and systematic (red squares) errors. The gray curve 
shows the fit with a p-wave Breit-Wigner function. 



tor accounts for the loss of available phase space for 
the production of the neutron and K + in the reaction 
p + p -> S(1385) + + K + + n with increasing A- 
7r + invariant mass and has been evaluated via simu- 
lations. The product of these two corrections is folded 
with the Breit-Wigner distribution lf26tl . The product of 
this Breit-Wigner distribution with the correction fac- 
tors is shown by the gray dashed curve in Fig. [7] The 
final result of the fit is also shown in Fig. [7] by the gray 
histogram representing the sum of all contributions. It 
delivers a normalized \ 2 value of about 1, a mass for 
the E(1385)+ resonance of 1383.2 ± 0.9 MeV/c 2 and 
a width of T = 40.2 ± 2.1 MeV/c 2 , where the er- 
rors are only statistical. Also other, commonly used, 
parametrisations of the Breit-Wigner have been tested 
and they have delivered results all compatible within 
statistical errors The signal, after the subtraction of the 
three background contributions, is shown in Fig.[8]to- 
gether with the relativistic p-wave Breit-Wigner func- 
tion ® obtained by fitting the distribution shown in 
Fig. |7] The experimental data are displayed together 
with the statistical and systematic errors. The system- 
atic errors have been evaluated modifying the analy- 
sis cuts as follows: the width of the cuts around the 



K + , A and neutron missing mass have been varied 
by ±20 %. Additionally, the cuts applied to recon- 
struct the A have been varied (d p _ 7T - < 28 mm and 
d(A — V) > 9 mm, see sectionUHl). For each set of cut 
values the same fitting procedure has been applied and 
a signal has been extracted after the background sub- 
traction. The obtained systematic uncertainty of the 
signal is shown in Fig.[8]by the red squares and lead 
to the following values for the pole mass and width 
of the resonance: m = 1383.2 ± 0.9±?;| McV/c 2 
andr = 40.2 ± 2.1 1*; 2 MeV/c 2 . The value of the 
width is about 4 MeV / c 2 larger than the PDG average 
quoted in ^ . On the other hand, the results obtained 
with this analysis take into account all the kinematic 
and efficiency effects that might modify the spectral 
function. The larger width of the S(1385) + resonance 
might also be due to its specific production mechanism 
in p+p collisions. Indeed, all the results collected in the 
PDG [3] refer to K~p and 7r + p reactions, where only 
a limited set of K~p data were used to arrive to the 
quoted average. 

Recently the CLAS collaboration published results 
on the £(1385) -> A + 7 H measured in photon- 
induced reactions; but no detailed analysis of the reso- 
nance mass distribution is reported upon. 



VI. DIFFERENTIAL CROSS SECTIONS 

The determination of the (differential) production 
cross sections requires extrapolation to full phase 
space. The corresponding (differential) acceptance 
corrections are obtained via a full scale simulation. Its 
input distributions are generated by a suitable model 
and tuned (iteratively) such that the model reproduces 
satisfactorily the obtained differential cross sections. 
The differential cross sections have been calculated in 
three reference frames which are described in the next 
section. 



A. Reference frames and angles 



The exclusive cross section in a 2 



3 reaction 



with unpolarized particles, do j \ C ^ L ^ L ^ L ), de- 
pends on four independent variables at a given value 
of y/s. In fact, for the reaction 1 + 2 — » a + b + c, 
the 9-dimensional exit phase space is constrained on 
a 5-dimensional hypersurface due to energy and mo- 
mentum conservation; azimuthal symmetry around the 
beam axis reduces the dimension of the exit phase 
space to four. Several choices of these four variables 
in different reference frames are conceivable 112811: 



• the opposite-momentum frame in the entrance 
channel (p[ — —pz) coinciding here with the 
1-2 center-of-mass system, 
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• an opposite momentum frame in the exit channel 
(Pa = -Pb)- 

In the first case, one may consider the angular distri- 
bution of the particle c with respect to the beam axis, 
yielding a &1_ 2 distribution. In the second case, one 
may consider the angular distribution of particle c with 
respect to the a — b direction yielding a Q c a _ b distri- 
bution or with respect to particle 1 yielding a 0*_ 6 
distribution. 0f_ 2 is usually labeled as CMS distribu- 
tion of particle c, Q c a _ b refers to the helicity angle and 
&a_ b is the Gottfried-Jackson angle. The CMS and 
the Gottfried-Jackson angles connect entrance and exit 
channels, while the helicity angles quantify relations 
within the exit channel. As the selection of particle la- 
bels a, b, c is arbitrary for different species, there are 
three different choices for each of the angles &1_ 2 , 
®a_ b and Q c a _ b - It has to be noticed that for identical 
particles 1 and 2, one should average 0*_ fe and Q 2 a _ b - 
For different hadrons a and b we supplement the super- 
script by the label of the hadron relative to which the 
angle is measured, e. g., Q^ b _ b or O^Z^ etc. 

The helicity angle distribution represents a special 
projection of a Dalitz plot. In particular a uniformly 
populated Dalitz plot results in an isotropic helicity 
angle distribution, whereas physical or kinematical ef- 
fects distorting the Dalitz plot reflect themselves in 
anisotropic helicity angle distributions. 

The motivation for an analysis within a Gottfried- 
Jackson frame arises from considering e. g. the distri- 
bution of the angle © E »Zjf-+> which is the angle of the 
£(1385)+ relative to the proton in the £(1385)+-# + 
reference frame. This angle can give insight into the 
scattering process especially concerning the involved 
partial waves. This statement holds true also if an in- 
termediate resonance A* is excited, e.g. np — > A* — > 
£(1385)+#+. 

B. Modeling of the acceptance correction 

a. Angular distributions in the center-of-mass sys- 
tem. To extract the differential cross sections as a 
function of the CMS angle, the data sample was di- 
vided into 7 angular bins and the analysis procedure 
described in sections [III] and [IV] was applied. This 
means that for each bin the neutron missing mass 
spectrum was evaluated and the background was de- 
termined, as described in section [IV] After the se- 
lection of the neutron, the £(1385) + signal was ex- 
tracted from the (A-7r + ) invariant mass distribution. 
The misidentification background was again fixed by 
the fit to the neutron spectra for each of the 7 angular 
bins. The background contributions due to the chan- 
nels p + p — > A + K + + n + 7r + and p + p — > 
S° + K + + n + 7r + were also fixed by using the scal- 
ing factor from the fit to the integrated spectrum (see 
Fig. [7]). In this way the background contribution could 
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FIG. 9: Color online. Angular differential cross section for 
the E(1385) + production as a function of cos(Qq MS ). The 
curve corresponds to a fit with a Legendre polynomials. 



be subtracted, and a pure S(1385) + spectrum was ob- 
tained for each bin. Applying this binning also to the 
full-scale simulations, the geometrical acceptance and 
efficiency could be extracted in a differential form and 
could be used to correct the differential distributions. 

The yield of S(1385) + in each angular bin was ob- 
tained by integrating the subtracted spectra in a ± 6a 
interval around the pole mass, where the parameters 
were taken from the fit to the integrated spectrum 
(Fig. [8j. The production cross section was normalized 
to elastic scattering which was independently mea- 
sured in the same experiment l32ll . Figure [9] exhibits 
the differential cross section for E(1385) + production 
as a function of the CMS scattering angle. The exper- 
imental data points shown in Fig. |9] point to a strong 
anisotropy for the emission of the S(1385) + . This 
anisotropy is quantified by fitting a series of Legendre 
polynomials with the coefficients A^, i= 0, 2 and 4 

— ^— = A Q ■ L Q + A 2 ■ L 2 +A A ■ L 4 , (4) 
a cost) 

Only even polynomials are used due to the symme- 
try in the entrance channel. The results of the fit are 
listed in the legend of Fig. [9] The non-zero contri- 
bution of the p-wave polynomials L 2 already reflects 
the peripheral character of the production mechanism. 
This differential distribution was obtained with accep- 
tance corrections determined by simulations in which 
particles are emitted according to phase space. This 
evaluation of the acceptance will not be sufficient if 
other, independent angular distributions show a devia- 
tion from an isotropic phase-space emission. 

The angular distributions of the neutron and K + 
were determined in the same way as for the §(1385). 
However, the interdependence of the three precludes 
further information about deviations from emission ac- 
cording to phase space. 
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b. Angular distributions in the helicity angle 
frame. We consider now the helicity distributions 
with respect to E _^-' tne an gl e between the 
£(1385) + and the neutron in the neutron-!^ refer- 
ence frame. For the acceptance corrections, phase 
space simulations filtered by the CMS angular distri- 
bution of the £(1385) + (see Fig. [9]) have been used. 
Since there is no obvious kinematical correlation be- 
tween the two angles 0^ and Q^MS' tne filtering 
is not expected to bias the extracted distribution in the 
helicity frame, but accounts only for the correct accep- 
tance of each event. As can be seen in Fig. [10] the 
helicity angular distribution is not isotropic but shows 
a peak around cos(8 *^_) = —0.5. This effect must 
be taken into account to get the appropriate acceptance 
corrections. 

Higher resonances, contributing to the production of 
the £[1385) + , could induce such an effect. In 
ref. 11411 . the production of £(1385) + for the same 
reaction system as analyzed in the present work, 
was investigated at a beam momentum of p = 
6GeV/c. They found that a part of the £(1385)+ 
production proceeds via an intermediate A ++ res- 
onance with a Breit-Wigner mass peak value of 
around 2035 MeV/c 2 , followed by the decay A++ -> 
£(1385) + + K + . The statistics was not sufficient to 
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FIG. 10: Color online. Angular differential cross section for 
the E(1385) + production as a function of cos(0^ A ,^ n ). The 
red and blue histograms correspond to the contributions by 
the non-resonant production of the E(1385) + and resonant 
production via an intermediate A ++ , respectively. 



obtain any information about the quantum numbers of 
this A ++ state, however, the width of the resonance 
was estimated to be about 250 MeV. These parame- 
ters were used as an input for a Monte-Carlo simula- 
tion of the process p + p — > A ++ + n with subse- 
quent A ++ decay into £(1385) + K + . From the sim- 
ulation we obtain the shape of the angular distribution 
of the £(1385) in the helicity angle frame. Figure [TOl 



shows, together with the experimental data, simulation 
results stemming from non-resonant production of the 
£(1385) + (red histogram) and via a A ++ production 
(blue histogram). The dashed curve shows the sum of 
the two simulated processes assuming a contribution 
of 66 % by the non-resonant production and 33 % by 
the A ++ excitation. One can see that the agreement 
with the data is excellent. This result suggests that a 
rather large amount of the extracted £(1385) + may 
be produced via an intermediate A ++ resonance. In 
the following it is assumed that indeed 33 % of the 
reconstructed £(1385) + stems from an intermediate 
A ++ . To obtain the appropriate filter functions for 
the simulations an iterative procedure is applied. In a 
first step the differential distribution of the £(1385) + 
in the CMS is used to filter the non-resonant simula- 
tions. The differential cross section of the neutron in 
the CMS is used to filter the A ++ simulations. This is 
a natural choice, as the A ++ and the neutron are going 
back to back in this reference frame. The acceptance 
corrections are recalculated using the filtered simula- 
tion and new experimental distributions are obtained. 
The process is repeated again and the new filter func- 
tion is extracted from the experimental distributions. 

The final results after three steps of iteration are 
shown in Fig. [TT] panel (a) for the £(1385) + distribu- 
tion in the CMS and in panel (b) for the neutron distri- 
bution in the CMS. In addition to the experimental data 
the two components of the simulated yields are drawn. 
The red curves represent the non-resonant contribution 
and the blue curve the A ++ contribution. The experi- 
mental data are shown together with the statistical and 
systematic errors (gray rectangles). The systematic er- 
rors have been evaluated by varying the analysis cuts 
in the following way. The cuts on the neutron missing 
mass and on the A invariant mass have been expanded 
by 20 %, the cut on the K + mass has been varied by 
+20 % and -10 %. The energy loss cuts in the MDC 
and TOFino detector have been also varied by ±10 %. 

The agreement of the simulated data with the ex- 
perimental points in panels (a) and (b) is a constraint 
imposed by the iterative procedure described above, 
on the other hand the validity of the new acceptance 
correction is supported by the fact that the distribution 
in panel (a) is more symmetric then the one shown in 
Fig. [9] In order to check whether the simulated model 
really fits to the experimental data, other, independent 
angular distributions have to be analyzed. Panel (c) 
of Fig. [TT] shows the CMS angular distribution for 
the K + which is also well reproduced by the simu- 
lated model. Panels (d), (e) and (f) depict the distri- 
butions of the three helicity angles and panels (g), (h) 
and (i) the distributions of the Gottfried-Jackson an- 
gles. Panel d) shows the distribution of cos(© s _^)_) 
already discussed in Fig. QJJ The small differences 
with respect to the previous figure arise by the mod- 
ified acceptance correction that now takes into account 
the CMS angular distributions. Panel (e) represents 
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FIG. 11: Angular differential cross sections for the E(1385) + production in CMS (top row: a : 0cmSi b : Qcmsi c '■ 

®cms)' helicity (middle row d : 6 S J^_, e : Q^tZ% + , f ■ 6j*^ + ) and Gottfried- Jackson angles (bottom row: g : 

n ~^L + , h : @v,Z-n) i '■ 0^»~^+) angle frames. The solid curves represent the contributions by the non-resonant (blue) 
and resonant (red) production mechanisms, as in Fig.[l0] The dashed curve shows the sum of the resonant and non-resonant 
contributions. 



the cos(02»_^ r+ ) distribution. There the anisotropy 
of the experimental distribution due the A ++ contri- 
bution is less evident but still present. The helicity 
angular distribution (panel (f)) is rather fiat for both 
the experimental and the simulated curves. A possi- 
ble anisotropy in the experimental data could be con- 
nected to the polarization of the A ++ that was not in- 
cluded in the simulations. The data leaves some room 
for anisotropy but no firm conclusion can be drawn. 

c. Angular distributions in the Gottfried-Jackson 
frame. The experimental distributions of the 



Gottfried-Jackson angles show a strong anisotropy 
that can be well reproduced by our simple model. The 
distributions shown in panels (g) and (i) of Fig. QT| 
are linked to the partial waves acting at the respective 
verteces but no specific assumption is made in our 
simulation. Nevertheless, these distributions are kine- 
matically correlated with the CMS distributions. This 
is also at the origin of the good agreement between the 
experimental and the simulated angular distributions 
in the helicity frame. 

It is clear that the nine angular distributions are not 
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all kinematically independent from each other, they are 
here shown for the sake of completeness. These distri- 
butions constitute an important reference to test future 
theoretical models describing the production. Indeed 
different production mechanisms will result in differ- 
ent angular distributions of the decay products. It has 
also to be pointed out that the internal degrees of free- 
dom of the S(1385) + and the intermediate A ++ reso- 
nance are neglected in our simulations. 

The overall agreement of the experimental data with 
the angular distributions modeled with the simulations 
justifies the usage of the resulting acceptance correc- 
tions, despite of the fact that the contribution by the 
A ++ resonance cannot be demonstrated unambigu- 
ously. 
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C. Production cross sections 

The last step of the analysis consists of the calcula- 
tion of the production cross section. The production 
cross section has been estimated by the integration of 
the simulated differential cross sections adapted to the 
experimental data over the whole phase space for each 
of the nine distributions shown in Fig. [TT] Indeed, as- 
suming that proper acceptance and efficiency correc- 
tions have been applied, the same total cross section for 
the exclusive production of the E(1385) + resonance in 
the reaction pp — >• n + K + + £+(1385) should be ob- 
tained. 

The integration of the nine differential angular dis- 
tributions shown in Fig. QT| delivers values of the total 
production cross section that are in agreement within 
the statistical errors. The arithmetic mean of these val- 
ues yields 22.42 ± 0.99 ± 1.57^^3 /zb. The statisti- 
cal error is followed by a first systematic error arising 
from the normalization to the elastic events and a sec- 
ond asymmetric error stemming from the cut variations 
discussed above. This cross section can be compared 
to the values known for the S(1192) + production in 
p+p collisions. Figure [T2l shows the production cross 
sections for the reaction p + p — >• E + + K + + n as 
a function of the excess energy e. The cross section 
extracted from our analysis of the channel p + p — > 
£(1385) + + K + + n is shown together with the mea- 
surement at higher energies reported in lfl2ll . 

The data point corresponding to the £(1385) + is 
about a factor two lower than the cross section value 
extracted for the E + at the same excess energy but 
still compatible within the systematic errors. Our data 
point agrees well with the value of 15 ± 2 /jb reported 
in (H for E(1385)+ at an excess energy of 830 MeV 
also in p+p collisions, pointing to a weak dependence 
on excess energy of both S + and E(1385) + produc- 
tion above e = 200MeV. 



FIG. 12: Color online. Production cross sections for the re- 
action p + p — > p + K + + E + as a function of the excess 
energy e. Data points measured by several experiments have 
been compiled in 1 33] . The production cross section for the 
reaction p + p — ► n + K + + E(1385) + from this work (red 
star) and at 6 GeV/c incident momentum 1 12] (blue star) are 
shown as well. 



VII. SUMMARY 

We present the results obtained from an exclusive 
analysis of the S(1385) + resonance produced in p+p 
collisions at a kinetic energy of 3.5 GeV. The sophis- 
ticated method developed for the reconstruction of the 
background allows to estimate with high precision the 
position and the width of this resonance. The value ex- 
tracted for the mass, m = 1383.2 ± 0.9^5 MeV/c 2 , 
agrees well with the PDG value 10], while the width, 
T = 40.2 ± 2.1±^MeV/c 2 , is about 4MeV/c 2 
larger than the PDG average yfl . 

Angular distributions have been corrected for accep- 
tance and detector response by means of a complete 
simulations testing two assumptions for the produc- 
tion mechanism. Our analysis suggests that a 33 % 
of the S(1385) + yield originates from the decay of 
an intermediate A ++ resonance, whereas one should 
underline that interference effects could also modify 
the kinematics of the reaction. The efficiency and 
acceptance corrected angular distribution of the three 
reference frames CM, helicity and Gottfried- Jackson 
have been extracted and can be used to test theoreti- 
cal models describing the production mechanisms for 
the S(1385) + state. A total production cross section 
of a = 22.42 ± 0.99 ± 1.57±|;^ £tb has been de- 
duced and is found to be consistent with the systemat- 
ics measured for the ground state S + production in the 
same final state as a function of the excess energy. The 
available data base do not exhibit a strong dependence 
of the production cross section on the excess energy in 
the range 200 < e < 800 MeV for both the ground 
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state S + and the first excited state E(1385) + . Our re- 
sults provide the necessary reference for further studies 
of the spectral shape of the £(1385) resonance in p+A 
and A+A reactions. Furthermore, our analysis repre- 
sents an important bench mark for the investigation of 
the A(1405) resonance in p+p, p+A and A+A colli- 
sions. Our results will help to understand the proper- 
ties and production characteristics of the lowest lying 
hyperon resonances both in elementary and heavy ion 
collisions with proton induced reactions at nuclei as 
link between them. 
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